Picosecond charge variation of quantum dots under pulsed excitation 
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We present a spectroscopic study of excitation dynamics in self assembled CdTe/ZnTe quantum 
dots. Insight into details of kinetics is obtained from the time resolved micro-photoluminescence, 
single photon correlation and subpicosecond excitation correlation measurements done on single 
quantum dots. It is shown that the pulsed excitation in energy above the energy gap of the barrier 
material results in separate capture of electrons and holes. The capture of carriers of different charge 
take place at different delay from excitation. 

PACS numbers: 78.55.Et, 78.67.Hc 



I. INTRODUCTION 

Quantum dots (QDs) belong to the most intensely 
studied topics in the solid state physics. They owe their 
popularity to new physics involved and to a wide range 
of their possible applications in such fields as fabrication 
of efficient light sources, single photon emitters, informa- 
tion storage and processing^ 2 -. A particular interest is 
related to the emerging field of quantum information 3 '^. 

Self-assembled semiconductor quantum dots receive an 
important share of the research effort, due to efficient fab- 
rication methods by modern epitaxial growth techniques 
and possibilities of integration with existing electronics. 
QD studies started from the prototypical InAs/GaAs ma- 
terial system, and were quickly extended over the entire 
families of III-V and II- VI semiconductors^* 6 -. Among the 
physical phenomena studied in the quantum dot research, 
those related to light emission represent an important 
part, both under resonant and non-resonant excitation. 

All-spectroscopic methods are well suited to study ex- 
citation and light emission processes in the QDs, in par- 
ticular their dynamics^. The most precise information 
on the physical mechanisms involved in the excitation 
and light emission processes is usually supplied by single 
QD spectroscopy. 

Basic QD spectroscopy methods include photolumines- 
cence (both cw and time-resolved) under varied experi- 
mental conditions: excitation power, temperature, etc. 
In time-resolved studies the temporal resolution is usu- 
ally determined by the type of detectors used (down to 
tens of ps for avalanche photo-diodes, to several picosec- 
onds for streak cameras). 

In some cases, more sophisticated spectroscopic tech- 
niques are necessary. For example, photon correlation 
measurements have been used to establish that in the 
case of non-resonant excitation of QDs, carriers are 
trapped separately rather than as whole excitons^ (sep- 
arate carrier capture in QDs was also demonstrated 
indirectly in cw experiment a 1 ^ 11 ). Standard form of 
pump-probe techniques is not much used, as absorp- 
tion measurements of QDs present serious experimental 
difficulties^^. 



A technique more feasible for QD studies, similar to 
but not identical with pump-probe methods, is excitation 
correlation spectroscopy (ECS)±i±£«±&±Z»i&±2. In ECS, 
photoluminescence is excited by pairs of laser pulses sep- 
arated by a controlled delay. It has been shown to be a 
powerful tool to investigate transient processes in semi- 
conductors, especially excitonic recombinatio n 16 ! 20 . It 
has advantage of outstanding temporal resolution lim- 
ited only by the properties of light pulses. However, not 
all the possibilities offered by ECS have been exploited 
so far. For instance, the order of carrier trapping in the 
excitation processes has not been studied to the best of 
our knowledge. 

In this work, we profit from the excellent temporal 
resolution of the excitation correlation spectroscopy and 
apply it to a study of population dynamics in single 
CdTe/ZnTe quantum dots (QDs). In particular we study 
dynamics of carrier trapping by a QD under non-resonant 
pulsed excitation. The choice of the CdTe/ZnTe system 
is motivated by its two advantages with respect to the 
classical InAs/GaAs one. First, light emission in the vis- 
ible range (red to green), and second, more robust ex- 
citonic states assuring efficient light emission at higher 
temperatures. 



II. SAMPLES AND EXPERIMENTAL SETUP 

The studied sample contained an MBE-grown single 
layer of self-assembled CdTe/ZnTe QDs. The sample 
growth was described in detail in Ref. [2l|. The density of 
quantum dots was estimated to be about 5 x 10 9 cm~ 2 . 
Measurements were performed on a sample immersed in 
superfluid helium (at 1.8K). A reflection microscope, im- 
mersed together with the sample, assured a spatial res- 
olution better than 0.5 /z:m. A frequency-doubled Sap- 
phire:^ femtosecond laser was used for pulsed above- 
barrier excitation. In excitation correlation experiments, 
the sample was excited by pairs of pulses with a con- 
trolled temporal separation (delay) between the pulses 
in a pair. Consecutive pairs were separated by the laser 
repetition period 13.6 ns. Time-integrated photolumines- 
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FIG. 1: Experimental setup used in excitation correlation 
experiments. The electro-optic modulator (EOM) was used 
to stabilize beam intensity probed after passing by a single 
mode fiber (SMF). BS denotes 50/50 beam splitter and P 
linear polarizer. 
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cence (PL) spectra were then recorded by a CCD camera 
as a function of the delay. The experimental setup is 
presented on Fig. [TJ 

Laser pulses were split in pairs in a Michclson inter- 
ferometer setup. The length of one arm of the inter- 
ferometer was varied by moving a corner-cube retrore- 
flector mounted on a motorized translation stage. The 
setup allowed us to achieve a controlled delay up to 4 
ns. Beams from two arms of the interferometer were 
combined again on a 50/50 beamsplitter forming a train 
of pairs of pulses. The joint beam was then transmit- 
ted through 50 cm of single-mode optical fiber acting as 
a spatial filter to assure a precise overlap of two laser 
spots on the sample. The width of each laser pulse at 
this point was estimated as 0.5 ps. The most challeng- 
ing task in such experiment was to assure the stability 
of the excitation of a single quantum dot. Due to imper- 
fections in optical alignment, the variation of the delay 
between the pulses led to changes in the efficiency of cou- 
pling to the fiber. This effect was canceled by introduc- 
ing an electro-optical modulator into the variable-length 
arm of the interferometer, to stabilize the intensity of the 
laser after the fiber. As a result, a good stability of exci- 
tation of a single quantum dot was maintained over the 
measurement time which could exceed 6 hours. 

In case of the single photon correlation measurements, 
a Hanbury-Brown&Twiss detection scheme^ was used. 
The photoluminescence from the sample was split on a 
50/50 beamsplitter and resolved by two monochroma- 
tors equipped with avalanche photodiode (APD from 
Per kin Elmer or IdQuantiquc) single photon detectors. 
The APDs were connected to START/STOP inputs of 
TimeHarp 200 time counting system. An electrical de- 
lay introduced in the STOP signal allowed us to detect 
photons at negative delay values. 



FIG. 2: (color online) (a) PL spectra of a single QD under 
different excitation intensities, (b) PL spectra excited by pair 
of laser pulses with temporal separation At. The At = Ops 
spectrum was shifted horizontally for clarity by —0.5 meV. 



III. PHOTOLUMINESCENCE SPECTRUM OF 
A SINGLE QD 

Microphotoluminescence spectra of QD ensembles, 
limited by the size of the excitation and detection 
spots, revealed an inhomogeneously broadened distribu- 
tion with a characteristic line structure. A low density 
of the lines in the low energy tail of the PL band allowed 
us to find well isolated sets of lines originating from sin- 
gle quantum dots. An example PL spectrum of a single 
quantum dot is presented on Figure O The lines were 
identified as originating from recombination of neutral 
and charged excitons and biexcitons, as marked in the 
figure. The identification was based on relative emission 
energies, in-plane anisotropy effects, and photon correla- 
tion measurements. 

The lines emitted by the QD in the neutral or singly 
charged state were first tentatively identified on the basis 
of the characteristic pattern of their emission energies, 
observed in previous experiments on similar samples 23 , 
in particular with charge-tuning 2 ^. The identification 
of the lines related to neutral exciton (X) and biexciton 
(XX) transitions was confirmed by characteristic in-plane 
anisotropy effects. In case of a quantum dot of C2v sym- 
metry, both X and XX lines are split in linearly polar- 
ized doublets, originating from the fine structure splitting 
(FSS) of the excitonic stated. The experimental resolu- 
tion did not allow us to observe the splitting directly. 
However, the components of the doublets could be ob- 
served in linear polarizations parallel and perpendicular 
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FIG. 3: (color online) Measured transition energy as a func- 
tion of detected light polarization. Oscillatory behaviour in- 
dicates a small anisotropic splitting. 



to the QD anisotropy axis. At intermediate polarization 
angles, each doublet was observed as a single broadened 
line at an intermediate spectral position. This effect leads 
to oscillations of the apparent X and XX transition en- 
ergies as a function of orientation of detection polariza- 
tion, as presented in Fig. [3] As expected, no energy 
oscillations were observed in case of charged excitons, in 
particular X + and X - , which contain pairs of identical 
carriers in singlet states. An argument supporting the 
assignment of trion signs is a negative optical orientation 
of X~ line, observed at quasi-resonant excitation through 
a neighbor do1i22. The negative optical orientation had 
been observed for negative trions in many QD systems 2 ^ 
and is related to electron-hole flip-flop process. 

Some linear polarization effects were also observed for 
the doubly charged exciton X 2- line. However, they are 
beyond the scope of this work^l. 

Further data supporting the identification of the lines 
were obtained from photon correlation measurements. In 
our experiments, the quantum dot was excited by pi- 
cosecond pulses of light and photoluminescence photons 
related to selected excitonic lines were counted by two 
detectors. Example results of such experiments are pre- 
sented in Fig. [U in the form of histograms of detection 
events of pairs of photons from the two transitions as 
a function of their temporal separation (number N of 
laser repetition periods). Due to pulsed excitation, time 
delay between the two emitted photons is close to integer 
multiples of the repetition period. A clear antibunching 
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FIG. 4: (color online) Typical photon correlations related to 
change of charge state. Negative time distance denotes pairs 
with inverted photon order, e.g. X - after X in case of (a). 
Solid lines were calculated within a model adapted from Ref. 
@ with parameters a = 0.80, f3 = 0.86, and £ = 0.26. 



(suppression of the peak) at zero delay confirms unequiv- 
ocally that we deal with a single photon emitter. Sim- 
ilar antibunching was observed for autocorrelation ex- 
periments, whereas cascade emission was witnessed by 
a characteristic bunching (enhancement of the central 
peak) in XX-X cross-correlation histograms (not shown) . 
The X 2 ~ line was identified using the cross-correlation 
histograms presented in Fig. Q] Besides the central anti- 
bunching, they show longer time-scale effects, extending 
over several repetition periods. Such effects are known to 
originate from QD charge variation 9 . The relatively high 
probability of the observation of correlated photons emit- 
ted after adjacent pulses in different charge states indi- 
cates an effective capture of single carriers. In particular, 
recombination of a neutral exciton after recombination of 
a negative trion requires only a single hole capture while 
a capture of three carriers is necessary if the emission 
order is opposite. Therefore corresponding probabilities 
of the events are respectively higher (N = +1 peak at 
Fig. QJi) and lower (JV = — 1 peak at Fig. (3^) than the 
probability in stationary state (N — ► oo). The similarity 
of X~ - X 2 ~ and X - X~ correlation histograms supports 
the assignment of the line X 2 ~ to the doubly charged 
exciton. This assignment is confirmed by relatively long 
characteristic time-scales of the X + - X 2 ~ correlation his- 
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FIG. 5: (color online) (a) Photoluminescence intensity of vari- 
ous excitonic lines and (b) estimated average QD charge state 
versus excitation (pulsed laser) mean power. Solid line was 
calculated within a model described in Section lV Bl and shifted 
vertically for clarity. 



togram, as it is related to the largest change of the QD 
charge. 

These qualitatively discussed correlation histograms 
can be simulated by a rate-equation model described in 
Ref. [^. We adapted this model by extending the possi- 
ble states to incorporate transitions in +1 and —2 charge 
states. Free parameters of the model include probabil- 
ities of capturing an electron, a hole or a free exciton. 
Solid lines in Fig. [5] are calculated with probability val- 
ues taken from Ref. 0. The results of the simulation 
confirm the predominant role of the single carrier trap- 
ping. They depend weakly on exact value of free exciton 
capture probability. A similar model is used in Section 
IV CI to simulate time profile of excitation correlation re- 
sults. 

The simplest experiment which gives an insight into 
the excitation dynamics is the measurement of photolu- 
minescence spectra at various intensities of the pulsed 
excitation. A typical dependence of the intensity of se- 
lected excitonic lines on excitation power is presented in 
Fig. [UJa). In the lowest excitation limit, intensities of the 
PL lines exhibit power-like dependence. In case of X, X + , 
and X~ transitions, the dependence is linear. Two inten- 
sities (X 2 - and XX) increase super- linearly. The superlin- 
ear dependence observed for the X 2_ line is not surprising 
in view of a large number (four) of carriers necessary to 
form the X 2 ~ complex. Most of these results are in agree- 
ment with typical behavior expected for single exciton 
complexes and for biexcitons when the capture of whole 
excitons is significant* 9 - ' 28 ' 29 . A further increase of the ex- 
citation power leads to saturation of the line intensities. 
This behavior is related to the fact that each excitation 



pulse results at most in one recombination cascade (the 
excitation time is much shorter than the recombination 
time). Therefore after one pump pulse only one pho- 
ton related to a determined transition may be emitted. 
Nearly quadratic dependence of biexciton PL intensity 
indicates a need for including process of free exciton cap- 
ture in addition to single carrier capture evidenced by 
photon correlation experiments. One should note that 
even a relatively small exciton capture rate may dom- 
inate over single carrier capture rates at low excitation 
power in multi-step excitation process, e.g. in case of XX 
formation. 

At medium and high excitation power, the lines corre- 
sponding to negatively charged states of the QD become 
relatively more intense. This effect can be analysed quan- 
titatively in terms of the average charge of emitting QD 
states calculated as: 



• 7 X + 1 • I x+ - 1 ■ J x - - 2 • I X 2- 



/x + / x + + h 



I X 2- 



(1) 



where Is is PL intensity of line S. This formula approxi- 
mates the averaged charge state of the quantum dot be- 
tween excitation events since only the fundamental tran- 
sition of each observed charge state is taken into account. 
No matter which was the highest state in a recombination 
cascade, its last step must be one of the final transitions: 
X, X + , X~or X 2 ~. Fig. [5]presents the average QD charge 
state as a function of excitation power. The effect of QD 
becoming negatively charged under strong excitation has 
been observed^ but the underlying mechanism cannot 
be determined without more detailed studies. It may be 
caused by a modification of the electrostatic environment 
of the quantum dot (similarly to that known for quan- 
tum well o 30 ' 31 ). However, a mechanism inherent to the 
quantum dot itself is also possible. It is related to the 
apparent absence of doubly positively charged states of 
the quantum dot, while doubly charged negative trions 
have been identified. The difference between number of 
bound states in positively and negatively charged QD is 
predominantly caused by the small valence band offset in 
the CdTe/ZnTe system. 



IV. MEASUREMENTS OF EXCITATION 
DYNAMICS 

Excitation dynamics of the QDs was studied by means 
of excitation correlation experiments. As explained in 
Section II, a selected QD was excited by pairs of laser 
pulses. Time-integrated PL intensity was measured as 
a function of temporal separation At between the two 
pulses in the pair. Plots of such dependence over the full 
temporal range are presented in Fig. [HIa). 

Features on two characteristic time-scales can be dis- 
tinguished: a dip several hundreds of ps wide, and a much 
sharper feature, both centered at zero delay. The main 
effect is the relatively wide dip in the PL signal. Its width 
is comparable to the radiative lifetime of excitonic states 
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(Fig. EJb)). The effect arises near the saturation regime 
when virtually each laser pulse excites the QD to a higher 
state. A qualitative explanation can be based on the fact 
that if the second pulse in a pair arrives prior to the 
excitonic recombination then the second pulse does not 
contribute to the intensity of the X transition. On the 
contrary, when pulses are separated by a few nanosec- 
onds, they act independently and therefore the recorded 
PL intensity is doubled with respect to single pulse exci- 
tation. Thus, in the simplest approach the dip should be 
described by an exponential function which may be di- 
rectly compared to the decay of the photolumincscence. 
This is shown on Figure [6{b-c) where the dashed line 
presents a profile obtained by fitting a monoexponential 
decay for |Ai| > 75 ps of the X profile from Fig. 6a. 

The second time-scale in the experiment is in the range 
of tens of picoseconds. An additional variation of the PL 
intensity is observed within this scale, as presented in 
Fig. \T[a). The signal is increased or decreased depend- 
ing on the excitonic complex with which the PL line is 
related. A clear increase of the photoluminescence at 
zero delay is seen for X~and X 2- . A decrease is seen for 
neutral exciton and X + . However, no significant effect 
is observed for a sum Epl of intensities of photolumines- 
cence lines related to all the observed charge states of the 
quantum dot (neutral and charged). This sum is shown 
on Figure [7] and its temporal variation is limited only to 
the slow component related to the excitonic decay. In- 
variability of the sum suggests that the effect is related 
to the QD charge state. Therefore wc analyzed the av- 
eraged QD charge state versus pulse separation. The 
calculated charge state is shown in Fig. [Tfb) and demon- 
strates a decrease of the averaged charge of the quantum 
dot when the pulses are in coincidence. The character- 
istic time of the feature appearing on the plots is the 
same as the lifetime of barrier luminescence, observed 
on similar samples 32 . Therefore we might associate this 
feature to the process of the carrier capture by the QD. 
The marked variation of the QD charge state during the 
pulse suggests a possibility of non-synchronous capture 
of carriers of different charge. In other words we wish to 
examine consequences of a delay between the arrival of 
holes and electrons in the QD. 

We parameterized the short-timescale feature using 
the following procedure. In the first step we sub- 
tracted a baseline originating from long-timescale effect. 
We assumed for simplicity the same characteristic long- 
timescale profile for each emission line and rendered it by 
a total PL signal Epl- We rescaled the total PL signal for 
each emission line by a constant factor to fit data points 
in a range 100 ps < |Af| < 125 ps. The rescaled tem- 
poral profile was then subtracted from the profile of the 
analyzed line. In the second step, we fitted an empirical 
function aexp(— \ At/b\) to the refined data, obtaining an 
amplitude a and a time constant b for each experimental 
scan. The results of this procedure are presented in Fig. 
EKc-d) versus excitation power. 
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FIG. 6: (color online) (a) Excitation dynamics over long time- 
scale. Presented data was symmetrized (averaged values for t 
and -t) for clarity, (b) Photoluminescence decay after a single 
laser pulse for neutral exciton recombination, (c) Effect of two 
pulse excitation in a semi-logarithmic scale. The dashed lines 
on all panels correspond to decay with 400 ps time constant. 



V. MODEL DESCRIPTION 

Experimental intensity-vs-delay plots, presented in the 
previous Section, exhibit two main characteristic fea- 
tures, centered around zero delay: (i) a sub-nanosecond 
decrease, common for most of the PL lines (with a char- 
acteristic time comparable to the PL decay time), (ii) 
a feature on the scale of tens of picoseconds, related to 
the variation of the QD charge state. A complete model 
description of the observed features is complex and re- 
quires certain assumptions concerning excitation mech- 
anisms of the quantum dot, its relaxation channels and 
their characteristic times. Therefore, to achieve a better 
insight in the physical mechanisms involved, we first dis- 
cuss simplified versions of the model, describing selected 
characteristic features of the data. In part (A) we discuss 
sub-nanosecond effects in the model neglecting details of 
the excitation process. In part (B) we show how a de- 
lay between the capture of electrons and holes results in 
a fast variation of the averaged charge of the quantum 
dot. The amplitude of this variation is described in a 
simplified model, in which both carrier capture profiles 
are completely separated in time, and only their integrals 
are meaningful. In part (C) we include an analysis of the 
temporal profiles and discuss characteristic times of the 
carrier capture. 



6 



CO - — - 
CO <2 




CD 

£ -1.0 

Q) 
> 
< 



-75 75 

Pulse 
separation (ps) 



0.0 0.5 1.0 1.5 2.0 2.5 
Excitation power (mW) 



FIG. 7: (color online) Picosecond scale photoluminescence 
dynamics: (a) example result of excitation correlation exper- 
iment for different excitonic transitions and their sum (b) av- 
erage charge state calculated according to eq. JT}. Power 
dependence of (c) amplitude and (d) time-scale obtained by 
fitting a |exp ( — At/6) | to baseline-corrected data (see text). 
Lines are calculated within a model described in Section V.B, 
with suitable x and y scale adjustments. 



A. Sub-nanosecond scale dynamics 
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FIG. 8: (a) Simulated temporal profiles for different excitation 
intensities fi. (b) The same simulated data in scale proper for 
a relaxation processes. Function exp (— t/400ps) is shown for 
reference (dashed line). 



The shape of observed long-scale PL dependence can 
be explained by introducing a simple analytical model. 
Here we neglect the effects related to the QD charge state 
and consider the QD energy spectrum as an infinite lad- 
der of states, starting from the lowest (ground) state for 
a given QD charge. Leaving out the charge degree of 
freedom in this model is justified by a similar decay dy- 
namics of all the observed states (Fig. d^b)). Within 
this simple model, we assume that the number of cap- 
tured e-h pairs after a single laser pulse does not depend 
on current QD state and is described by a probability 
distribution R(k). For example, in case of free exciton 
trapping the probability of capturing exactly k e-h pairs 
is given by a Poisson distributio n 28 ' 29 ! 33 

p(fc lA i) = ^-e^ and V(0,0) = 1 (2) 

where the average excitation fi is identified with the exci- 
tation power. In case of separate capture of electrons and 
holes in a dot with limited number of charge states, the 
probability distribution is close to the square of the Pois- 
son distribution. The particular shape of the distribution 
does not change dramatically the sub-nanosecond effects. 
After the pulse, the QD relaxes towards the ground state 
by radiative decay. We assume for simplicity that all 



the excited states have the same lifetime of r = 400 ps. 
Therefore the probability, that the quantum dot excited 
to n th excited level emitted exactly I photons over time 
t is given by truncated Poisson distribution: 

( V{l,t/r) if I < n 

■P n (l,t/r)={ EZnnht/r) if i = n (3) 
[ if / > n 

Within this model, we derived the following expression 
for PL intensity of the first excited state (i.e. exciton 
state) in an excitation correlation experiment, when the 
two pulses are separated by At: 

/(At) = (l-R (0)) ^1 + £ R (i) Vt (i, At/rU (4) 

Profiles of intensity versus At, simulated for different ex- 
citation intensities /i, taking R (k) = V(k,fi), are pre- 
sented in Fig. 8. 

It is interesting to note that the form of the simulated 
profiles is approximately mono-exponential (Fig. 8b). 
However, their slope decreases with increasing excitation 
intensity, and follows the lifetime of the system only in 
the limit of low excitation intensity. This finding should 
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be kept in mind when interpreting the measurements in 
which correlation excitation spectroscopy is used for de- 
termination of lifetimes 16,20 . The stability of our exper- 
imental setup was not sufficient to verify the predicted 
slope variation. On the other hand, the assumption of 
equal lifetimes of all contributing excited states is not 
fulfilled. 

In spite of the simplicity of the approach, in which we 
neglect any possibility of the change of the charge state of 
the quantum dot, the above model reproduces quite well 
sub-nanosecond temporal profiles obtained in the excita- 
tion correlation experiment. It gives also a quite good 
prediction of the variation of the X and XX intensities 
versus excitation power up to the saturation level. The 
predictions of the model for X and XX lines are shown on 
Figure 5a by solid and dashed lines respectively. The in- 
troduction of the possibility of the change of the charge 
state of the quantum dot requires taking into account 
details of the carrier capture after the excitation pulse. 



B. Picosecond scale dynamics 

The understanding of the short-scale PL dynamics re- 
quires a different approach. As mentioned previously, the 
effects on this scale are related to the QD excitation pro- 
cess (mainly the QD charge variation) rather than to the 
radiative relaxation after the first pulse. 

Generally speaking, the excitation correlation signal 
results from a non-additive character (non-linearity) of 
the QD excitation by the two pulses. This non-linearity 
may be attributed to processes occurring in the barrier or 
in the QD. Here we assume that both pulses generate the 
same number of carriers in the vicinity of the QD, and 
the non-linearity originates in the QD itself. To describe 
the observed QD charge state variation, we assume that 
electrons and holes exhibit different trapping dynamics 
after a single excitation pulse. We will show that this 
approach allows us to explain in a simple way all the 
observed experimental features. In particular, we repro- 
duce qualitatively the variation of the averaged charge 
of the quantum dot with excitation parameters: excita- 
tion power, delay between light pulses, and the ratio of 
energies of two consecutive light pulses. 

The non-synchronous trapping of holes and electrons 
may be caused by different processes which result in dif- 
ferent temporal profiles g e (t) and gh{t) of their capture 
rates. At this first stage, we will discuss a simple over- 
drawn case. We consider the carrier capture profiles as 
non-overlapping narrow pulses, with the electron trap- 
ping pulse delayed by time T e -h with respect to the trap- 
ping of holes. The carrier trapping pulses can be de- 
scribed by standard rate equations 



P+ = ~9eP+ + ghPo 

Po = 9eP+ - (gh + 9e)P0 + 9hP- 

P- = 9eP0 ~ (9h + 19e) P- + gh.Pl- 

P2- = 19eP- ~ 9hP2- 



(5) 



where we introduced an additional parameter 7 to ac- 
count for electron-electron blocking. The value of 7 was 
set at 0.3 on the basis of relative intensities of neutral- 
to-charged exciton line. Under our assumptions, the rate 
equations can be integrated separately for the electron 
and hole capture pulses, producing matrices A and B, 
describing the influence of the pulses on the charge state 
probabilities. 



= exp 
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■ 
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r 











-r 
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-r. 



= exp 



-r 
r 









-r 

r 








- 7 r 
7 r 

where T denotes the integral of capture rate over the 
pulse, assumed to be equal for electrons and holes and to 
be proportional to the intensity of the laser beam. 

Each pair of laser pulses in the excitation correla- 
tion experiment will produce two pairs of carrier cap- 
ture pulses. Depending on the separation between the 
two laser pulses, the hole trapping after the second pulse 
will take place before or after trapping of electrons from 
the first pulse. Thus, there are two possible orders 
of carrier trapping: holc-hole-electron-electron or hole- 
electron-hole-electron. The two cases can be described 
in terms of a recursive equation that binds charge state 
distributions before and after a pair of excitation pulses. 
For example, for the hole-electron-hole-electron ordering: 



(after) 
P+ 
(aftcr) 

Po 

(after) 
p_ 

(after) 

P2- 



A • B • A ■ 



(before) 

p\ 

(before) 

Po 

(before) 
(before) 

P2- 



(7) 



We computed stationary states in both cases and used 
the difference between them as a measure of the ampli- 
tude of the picosecond-scale feature. The results of the 
simulation are compared with the experimental results in 
Fig. [7]Jc), after an appropriate adjustment of both am- 
plitude and power scales. A good agreement is achieved 
at low excitation power, while some discrepancies appear 
at higher power. They may originate from the absence of 
higher charge states in the model description. The model 
provides also a correct qualitative description of the ob- 
served evolution of the average QD charge state towards 
more negative values under increasing excitation power, 
as shown in Fig. OJb). The enhancement of the negative 
QD states is obtained only if r e _^ > 0, that is if electrons 
are captured after holes. 

An additional test of the model was provided by exper- 
iments in which the intensities of two consecutive light 
pulses were different. This eliminates the temporal sym- 
metry of the obtained profiles. In our experiments, the 
power of one beam was kept constant and the power 
of the second one was set at different levels in a series 
of consecutive measurements. Typical experimental pro- 
files, obtained for the power ratio 1:2, are shown on Fig. 
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FIG. 9: (color online) (a,b) Results of excitation corre- 
lation experiment with two pulses of different intensities 
(0.5/iW/l.O/iW denotes pulse order for negative pulse sep- 
aration). (c,d) Plots of effect asymmetry obtained by sub- 
traction: PL(t) — PL(—t) for three sets of pulse intensities: 
0.2/0.5/nW, 0.5/0.5/iW, and 1.0/0.5/iW. 
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FIG. 10: (color online) Amplitude of asymmetry plotted 
against intensity of one of the pulses. Intensity of the other 
pulse was set to 0.5 /xW. Solid lines are calculated within 
a model described in Section V.B and rescaled by a common 
factor. 



[9ja). The obtained profiles are clearly asymmetric. This 
asymmetry is better visible in the plots of difference be- 
tween signals measured at opposite delays. An exam- 
ple is presented in Fig. [21(c) for different power ratio 
values. At ratio 1:1, as expected, the signal is almost 
zero for all the excitonic complexes. For ratio different 
than 1 , the lines related to negative exciton complexes are 
stronger when the stronger pulse arrives second. At the 
same time, neutral and positively charged exciton lines 
are less intense. The signal asymmetry increases during 
the first tens of picoseconds, reaches a maximum for a 
delay of about lOOps and then decays with a decay time 
similar to the exciton recombination time. We extracted 
the asymmetry amplitude and compared it to the predic- 
tions of the model. The amplitude obtained for different 
exciton complexes and ratios of the light power in two 
beams is presented on Fig. 1101 and marked by symbols. 
In the model, different powers of the two beams were 
simulated by taking different values of T (proportional to 
beam power) in two pairs of matrices A and B in equation 
7. The results of the simulation are presented by lines on 
Figure ITD1 The agreement is quite good and supports our 
interpretation of separate capture of electrons and holes. 



C. Continuous rate equation model 

The simple model discussed in section B was sufficient 
to analyze the amplitudes of the picosecond scale fea- 
tures observed in the excitation correlation experiment. 
However, to describe the temporal shape of the observed 
peaks, we need certain assumptions about the profiles 
of the carrier capture rates. Our data do not allow us 
to determine exact profiles, nevertheless they give some 
insight in the characteristic times. We propose here a 
rate equation model with simple exponential decays of 
the hole and electron capture rates gh(t) and g e (t). They 
both start at the time of arrival of the light pulse and de- 
cay with different time constant r/j and r e respectively. 
Such profiles could be related to the exponential decay of 
free carriers in the barrier material and/or wetting layer. 
The free carriers could be trapped by quantum dots or 
other centers. Direct measurements of the time resolved 
photoluminescence from the barriers in similar samples 
show a fast monoexponential decay^. Such decay would 
be a straightforward consequence of the above assump- 
tion if r e 3> Th, and then the measured PL decay was 
equal Th- 
in our continuous model, we consider 10 states of the 
quantum dot including ground and first excited states 
for the total QD charge +1, 0, -1, -2, and second ex- 
cited states for the charge and -1. This selection is 
based on the identification of the optical transitions ob- 
served in the spectrum. Following Ref. [^, we assume 
excitation by trapping of an electron, a hole or an entire 
exciton. Relative integrated rates of these processes are 
taken from Ref. QJ. Due to the double pulse excitation, 
the temporal profile of excitation of each type is a sum of 
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FIG. 11: (color online) Model of time dependence: (a) Cap- 
ture rates for both carrier types after a single laser pulse, 
(b) Calculated photoluminescence of X 2- and X + lines for 
various delays in two-pulse experiment, (c) X 2- simulations 
for r e = 25 ps and different parameters th- The curves are 
compared with experimental data (symbols) and simple ex- 
ponential relaxation with time-scale r e + 7>j. 



two exponential decays starting at arrival of subsequent 
laser pulses. The rate equations that include thus defined 
excitation and radiative recombination are integrated nu- 
merically. Photoluminescence intensities of various lines 
are found by integrating the respective radiative recom- 
bination over one repetition period, after finding a steady 
state of the system. 

Example temporal profiles, calculated for the extreme 
charge states (X + and X 2- ), are presented in Figure 
ITTT b). A sharp increase of the intensity of X 2 ~ is ac- 
companied by a decrease of the X + intensity. Both fea- 
tures have a similar shape which can be approximated by 



an exponential function. The characteristic time of this 
function is determined by parameters Th and r e . However 
it is never smaller than the larger of them and it is close 
to the value of their sum : Th + T e . The amplitude of 
the peak decreases when values Th and r e become closer 
to each other. A comparison of the characteristic times 
obtained from simulations to the experimental ones al- 
lows us to conclude that the electron capture time (the 
longer one) is in the range of 20-40ps. The hole capture 
time is much smaller to assure a sufficient amplitude of 
the observed features. 



VI. SUMMARY 

We performed a detailed time-resolved spectroscopic 
study of single CdTe/ZnTe quantum dots. The exci- 
tation dynamics was investigated by the time resolved 
micro-photoluminescence, single photon correlation and 
subpicosecond excitation correlation measurements. The 
time resolved experiments were done with excitation by 
ultrafast pulsed laser working at energy above the energy 
gap of the barrier material. The obtained temporal pro- 
files of excitation correlation exhibit several characteristic 
features: sub-nanosecond decrease of the intensity com- 
mon for most of the PL lines (with characteristic time 
comparable to PL decay time), and picosecond variation 
of the relative intensity of the lines related to excitons 
of different charge state. We propose a model describing 
observed features and demonstrate that it requires that 
the carriers are trapped separately. Moreover the cap- 
ture of carriers of different charge take place at different 
delay from excitation. The detailed analysis of the tem- 
poral profiles let us conclude that the electron capture 
takes place in 20-40ps after excitation, while capture of 
hole is much faster. 
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